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Abstract 

The Green's function for a hole moving in an antiferromagnet is derived ana- 
lytically in the long- wavelength limit. We find that the infrared divergence is 
eliminated in two and higher dimensions so that the quasiparticle weight is fi- 
nite. Our results also suggest that the hole motion is polaronic in nature with 
a bandwidth proportional to t/Jexp[— c(t/J)^] (c is a constant). The connec- 
tion of the long-wavelength approximation to the first-order approximation in 
the cumulant expansion is also clarified. 
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It has been recognised for some time that understanding the motion of a hole in a 2D 
antiferromagnet (AF) would be an important first step towards a complete understanding 
of the effect of doping on the Cu02 planes of high Tc cuprates, which are known to show 
antiferromagnetism in the undoped case. The AF with one hole is also a highly non-trivial 
correlated electron system, and is therefore of fundamental interest from a purely theoretical 
point of view. There have already been many studies of the one-hole problem including 
those based on exact diagonalizations (ED) of small clusters the self-consistent Born 

approximation [|6|-|T0[1, the restricted basis method the Bogoliubov-de Gennes (BdeG) 



equation ||T2[, and classical descriptions 0,0] (see |l| for further references) 



The few analytical results which do exist have proved very valuable. However, most of the 
previous studies have involved numerical calculations (even the studies using the SCBA have 
to solve Dyson's equation numerically for small clusters). This has left a few points which 
still need clarification. We mention two of these. Firstly, in the small J/t limit the SCBA 
yielded a power law dependence for the hole band width [~ t(J/t)°] |^. This is consistent 
with the results of ED's on small clusters and gives support to the 'string' picture [|T|. 
Numerical calculations based on the BdeG equation and another variational approach 
^ have suggested that the motion of the hole is polaronic in nature in a wide parameter 
region. However, both of the approaches use a Born-Oppenheimer type approximation, 
which explicitly breaks translational invariance, and it is not clear how much this may have 
affected the conclusions. Secondly, although numerical calculations on clusters show that 
the hole has a finite quasiparticle weight, there is still some uncertainty as to whether the 
quasiparticle weight vanishes or not in the thermodynamic limit [|15] . 



Here, by treating spin waves in the long- wavelength (continuum) limit, we derive an 
analytical expression for the hole Green's function for the hole momentum close to the band 
minimum (7r/2,7r/2). As we work directly in momentum space, translational invariance is 
not broken. This allows us to confirm explicitly the polaronic behaviour of the hole. The 
analytical expression we obtain for the hole Green's function can also be used to directly 
examine the hole quasiparticle weight. Our expression for the Green's function shows that 



the infrared catastrophe, which leads to the vanishing of the quasi-particle weight in the 
ID case, is ehminated in 2D (and higher dimensions), so that there is a finite quasiparticle 
weight in 2D. This is consistent with many other studies. We also show that our approach 
is equivalent to a cumulant expansion and is suitable for large and intermediate J/t, as for 
the usual polaronic problem. Our approach is therefore complementary to the SCBA, which 
is better in the small J/t limit. 

Our study is based on the t—J model. Treating the spin waves as the collective excitations 
of the Heisenberg antiferromagnet, the following effective Hamiltonian for the t — J model 
has been obtained by previous authors 

Hi = Ho + V, 

Ho = Y.^,P% (1) 

tz 



Ntg 



Here and I3g are the annihilation operators of the hole and the spin wave, z is the 
coordination number (z = 4 for a 2D square lattice), 7^ = Yl,h^^^'^ 1^ with 5 the unit vectors 
to nearest neighbors, and is the number of the lattice sites. The spin wave excitation 



spectrum ujg = JzsVg with Vg = y^l — 7^ and s = 1/2. The Bogoliubov transformation 
coefficients are Ug = [(1 + z/q)/(2z/g)]^/^ and Vg = — sgn(7g)(M^ — 1)^/^. Although (|l|) is not an 
exact mapping of the t — J model, very good agreement between the results obtained from 
the effective Hamiltonian and those from the original t — J model have been demonstrated 



for small clusters by many authors |]8|-|T0|. We take the Hamiltonian (1) as our starting 
point. 

The hole Green's function is defined as: 

G{k,i) = -z{Th^{t)hk\0)), (2) 

where t denotes time throughout the paper, h^it) is the Heisenberg operator with respect 
to H. The thermal average (■ ■ ■) is for the spin subsystem. Since there is no hole for the 
spin subsystem, we can write 



oo 

= -t9ii)Y.Gm{k,i), (3) 



m=0 

where 

^2m ri 



Gm{k, t) = dti ...j^ dt2m{Thkii)Vih) ■ ■ ■ VihmH). 



(2m)!iV^. ^ 



fdh . . . f di2n.{TMk,^,g,Jhm) ■ ■ ■ Mk,,,Sl)) 
u „ -10 Jo 



Here Pk,q(t) = h\_^{t)hk{i), and 

Mk,q{t) = {Uq-ik-q + Vq-fk)Pl{i) + {Ug^fk + Vg-fk-q) P-q{i) ■ (5) 

Formally we are treating as a perturbation. The operators 0{t) are now defined in an 
interaction picture with: 0{t) = exp{iHot)0 exp{—iHoi) . 

The hole part in (|), {0\hk{t)pk,^^q,Jt2m) ■ ■ ■ Pk^,qAii)hl{0)\0), equals (0|/ifc(t)4(0)|0) 
when ki = k — Y4Z1 Qh and is zero otherwise. We can therefore trace out the hole part and 
write @ into 







/ dt2m. ■ ■ ■ 


/ dh 
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■■■Mk-q,,. 



A^™(2m) 

(6) 
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In general, it is impossible to obtain an analytical expression for for large m since the 
momentum ki in Mk^^q^iti) varies with time tj, refiecting the "history" of the distortion of 
the spin background induced by the hole. 

Here we derive an analytical expression for the Green's function treating spin waves in 
the long wavelength or continuum limit. Based on Wick's theorem, the expectation value 
(see (D), (Mfc2^^g2^(t2m) ■ ■ ■ Mk^qSi))i is composed of contractions like 

Bkikjq,q,(ti,ij) = {Mk^^qXii)Mk^^q^{ij)) 
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where 



(8) 



We find that 
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{qi^ sin + qiy sin kyf + 0(g^, cos^ k) 



(9) 



for /c near the band minimum /cq and for small momentum transfers q. The important feature 
of Eq (H) is that Ckikjqiq^ does not depend on the exact value of the hole momenta, ki and 
kj, but only on the momentum transfer and the initial momentum k. This result requires 
only that 7^ = 1 + O(g^) for small q and, for k = kQ + {dk^, Sky), 7^ = {Sk^ + 5ky) + 0{5k^). 
We can then formally rewrite Gm{k, t) in (D as 



Gm{ki i) 



(itz) 



{TMk,q,Jhm)---Mk,q^ih)) 



dt 



2m 



(10) 



where 



Mk,qit) 



Using d^), we have 



V^q 



IqxSink^ + qySmky\[Pg{i) +/?l,(t)]. 



(11) 



G{k,t) = -ie{t){Te 



) 



(12) 



and hence 



G{k,i) = -ie{i)e 



-iekt - (pkii) 



(13) 



with 



Cfc = - 



•Xt) = - 



1 / t \ 2 1 g:;. sin fc,. + g,, sin ^ , . r , , , , • f 

- (^-j - ^ ^[(e— ^* - l)(iV, + 1) + (e-^* - l)iV,]. 



(14) 
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Here we will only discuss properties of the hole at zero temperature. At T = 0, A^^^ ^ 0. 
In ID, after performing the sum (integration for an infinite system) over q in (p!^), we 
find that the function (l)kiJ)\k=-K/2 ~ lii(l + iS,t) here t = i/Jzs and the momentum cut-off 
^ ~ 0(1). We have approximated the spin spectrum by linear dispersion. The logarithmic 
divergence at large times leads to the so-called orthogonality catastrophe. The corresponding 
spectral function A{ko,u) ~ 9{u — eko){uj — e^o)^^^, with g = (t/J)^/27r, and there is no 
quasiparticle behavior in this case. In 2D, however, we obtain that (pkoii) ~ it~^(e~*^* — 1)— ^. 
For large i (or t), since the first term is irrelevant, there is no logarithmic term. The constant 
term left in (pko {t) at the large t limit contributes a finite quasiparticle weight and the spectral 
function A{k,uj) is then of the form 

A{k,uj) = 2nZkS{uj - e^) + Anc(^,t^), (15) 

where the quasiparticle weight = exp[—Ck(t / J)^] [ck are constants). This exponential 
factor is reminiscent of the Huang-Rhys factor |T7|JT8| in the usual electron-phonon problem, 
indicating a polaronic behavior. This Huang- Rhys factor is in agreement with that obtained 
using the BdeG equation. |T2| 



Using (|IB|) and linear dispersion as the spin wave spectrum, we obtain that in two 
dimensions the hole spectral function at ko is given by 

A{ko, uj) = 27tZ6{lj - efc) + Re H dU^^-'^^^~' Z U'^i^'^^^' I't _ ij ^ (15) 

JO L J 

where c = Jzs and a = (t/ J)^/(47rc). The second term in (|1^), i.e. the incoherent part, 
is well-behaved. A{k,uj) is shown in Fig. 1. The incoherent part is almost constant over 
a broad energy region. In the results of very small cluster calculations (both exact and 
SCBA) p|,|2T|, many secondary peaks in A{k,ij) were found above the lowest quasiparticle 
one. These secondary peaks were attributed to "string" resonances. However, the cluster 
calculations also show that, when the size of the system increases, these peaks become less 
prononced and, recently, Leung and Gooding 0] found in exact diagonalizations that the 
secondary peaks which are well defined for a 16-site lattice disappear in a 32-site system 



and that the secondary peaks are a finite size effect. Our results are consistent with this 
suggestion. If the secondary peaks were smeared out, the spectral function obtained from 
the small-cluster calculations would have the same broad feature as that shown in Fig. 1. 

We note that the hole Green's function (^) has the same form as that of the first-order 
approximation in the cumulant expansion. In fact, in the cumulant expansion IT? 



Gik,t) = -teit)exp{J2F^{k,t}}, (17) 

n=l 

where 

F2ik,i) = G2ik,i)-^Ff, (18) 
F3{k, t) = Gsik, t) - F,F2 - ^Fl 

Here Gi is defined in (|^). The first order term is given by 

Fi(fc,t) = {tzf ^^^^^-^ + {zt + [(e-^-^* - l)(iV, + 1) + (e^-«* - l)iVj /oo,} . (19) 

q 

Since in the long- wavelength limit the spin waves are uncorrelated for the hole at the band 
minimum, Fi{k,i) = for i > 2 and Fi{k,i) = —ieki— (pkif)- So the result of the first- 
order approximation in the cumulant expansion is exact in the long-wavelength limit for 
the hole momemtum k = ko. For the usual polaron problem |T9|j20[] , numerical calculations 
have indicated that the cumulant expansion converges rapidly for weak and intermediate 



couplings |2^. But quite why the resummation into the exponential like the cumulant 
expansion is a proper choice for the problem has not been understood clearly. Here for the 
spin polaron problem, we establish a connection between the first order approximation in 
the cumulant expansion and the long- wavelength approximation for the boson excitations. 

The LWA gives that the energy of the hole scales as t'^/J, as seen in (|1^, which should 
be correct only in the large J/t case. In the intermediate and strong coupling regions, the 
BdeG study suggests that the energy should be multiplied by the Huang-Rhys factor. 



in which case the bandwidth should be represented by 



W/t = a-e-'^''^^\ (20) 

where a and c are independent of t and J. To see how well this universal functional can 
represent the hole bandwidth, we compare it with numerical results of other studies. We 
show various estimates for the hole bandwidth in Fig. 2. The dashed-dotted line and 'V 



correspond to results obtained using the variational approach by Sachdev and ED on 
clusters of 20 sites by Poilblanc et al 0. The variational approach is reliable in the large 
J/t case, while most ED's are only available for small J/t. We choose a = 2.8 and c = 0.5 
for the functional dependence to fit results of other studies. The functional dependence is 
shown by the solid curve in Fig. 2. It is very close to the variational results for large J/t and 
to the ED results for small J/t, especially J/t > 0.4. For small J/t, a power law fit to the 
functional dependence gives b + d{J /t)'^'^^'^ (the dashed curve). The coefficients h and d are 
different from those obtained from the numerical calculations The form (^) is thus 

only qualitatively correct in the small J/t limit. In the region of J/t between 1.0 and 2.0, 
the functional dependence describes a smooth crossover from t/J behavior in the J/t limit 
to roughly J/t behavior in the small J/t limit. For comparison, the results of SCBA (open 
circles) are also shown in Fig. 1. The SCBA seems to substantially underestimate values of 
the bandwidth in the intermediate and the large J/t cases p. 

We would like to mention that although the energy of the hole in (|1^, which is 
proportional to sin^ + sin^ ky, gives rise to the hole pockets at k^ correctly, the other 
terms like cos k^ cos ky, which lead to anisotropy of the effective mass, are of higher order 
in LWA. So the LWA cannot be expected to give the anisotropy quantitatively. To discuss 
the anisotropy of the effective mass qualitatively, one could use the result of the first order 
approximation of the cumulant expansion, i.e., Eq. ([T9|) . 

In conclusion, we have derived an analytical expression for the Green's function of the 
hole moving in an antiferromagnet near the band mimimum in the long wavelength limit. 
The Green's function clearly indicates that the infrared divergence is eliminated in two di- 
mensions so that the quasiparticle weight is finite. It also suggests that the hole motion has 



8 



a polaronic nature for intermediate and large J/t. We have shown that the cumulant ex- 
pansion is a good choice for studying the hole motion in the weak and intermediate coupling 
cases, with the first-order approximation equivalent to the long-wavelength approximation 
at the band minimum of the hole. This should be complementary to the self-consistent Born 
approximation which is better for small J/t limit. 

One of the authors (YML) acknowledges support from the EPSRC of the United 
Kingdom under grant No. GRK42233 and from MURST/British Council under grant No. 
Rom/889/92/47. 
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FIGURES 

FIG. 1. The hole spectral function A{ko,uj) as a function of u>. We take ^ = vr and J/t = 1.0. 
uj is in unit of Jzs. 

FIG. 2. The hole bandwidth W/t as a function of J/t for a hole moving in a 2D antiferromag- 
netic background. The functional dependence W/t = 2.8 jc"^'^^*/"^^^ (see the text) is shown by the 
solid curve, along with the results from the exact diagonalization calculations on 20 sites ("*'') 
from the self-consistent Born approximation on a cluster of 16 x 16 (open circles), ^ and using a 
variational approach. [^2| The dashed curve, which is proportional to (J/t)^'^'^'^ , is the best fit to 
the solid curve for small J/t. 
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